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OBJECTIVE 

Develop  an  approximate  model  for  the  vibration  of  a  two-layered  plate.  This  model 
will  be  applied  to  the  theoretical  investigation  of  the  scattering  of  acoustic  waves  from  sub¬ 
merged  structures  with  coating  layer(s). 


RESULTS 

1.  A  single  equation  is  developed  for  the  transverse  motion  of  a  two-layered  plate. 

2.  An  equation  is  derived  for  the  phase  velocity  of  straight-crested  waves. 

3.  Phase  velocities  of  the  first  four  antisymmetric  modes  of  motion  are  calculated 
for  an  illustrative  case  of  a  steel  layer  and  a  copper  layer. 


RECOMMENDATIONS 

1.  Employ  this  model  to  investigate  the  behavior  of  the  scattered  signal  (target 
strength)  from  a  two-layered  plate.  A  useful  approach  is  to  assume  the  plate  to  be  loaded 
with  water  on  one  side  and  air  or  water  on  the  other  side,  then  to  compute  and  plot  the 
transmission  and  reflection  coefficients  of  an  incident  acoustic  wave.  The  procedure  for 
performing  the  calculations  is  very  similar  to  that  of  Graff  et  al.  The  variation  of  the  coeffi¬ 
cients  versus  the  angle  of  incidence  and  the  wavenumber-thickness  product  yields  informa¬ 
tion  about  target  strength  and  its  variation  at  various  modes  of  plate  vibration.  Furthermore, 
since  the  coefficients  can  be  computed  from  the  exact  theory,  a  quantitative  comparison  can 
be  made  and  the  validity  of  this  approximate  model  can  be  examined. 

2.  Apply  this  model  to  a  composite  plate  of  one  elastic  layer  (base  plate)  and  one 
viscoelastic  layer  (coating  layer),  to  calculate  target  strength.  Impedance  discontinuities 
along  the  plate  such  as  stiffening  ribs  and  varying  thickness  can  be  handled  by  appropriately 
modifying  the  model. 

3.  Note  that  this  model  can  also  be  used  where  the  viscoelastic  layer  contains 
many  inclusions  and/or  cavities,  provided  the  “equivalent  effective  elastic  properties”  of 
the  viscoelastic  layer  are  calculated  and  used. 

4.  If  experience  with  a  single-layered  plate  is  any  indicator,  this  model  should  be 
improved  to  provide  better  resolution  of  the  scattered  signal  of  acoustic  waves  at  grazing 
incident  angles,  possibly  by  combining  it  with  the  “Lyamshev  Theory”  on  symmetric  modes 
of  vibration. 
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INTRODUCTION 


When  an  acoustic  wave  strikes  a  plate  of  homogeneous,  isotropic,  and  elastic  material, 
the  plate  vibrates  in  two  manners  of  motion:  extensional  and  flexural,  each  containing  infi¬ 
nite  numbers  of  modes.  In  extensional  motion,  the  displacement  of  the  material  particle, 
averaged  over  the  thickness,  is  in  the  direction  parallel  to  the  plane  of  the  plate.  This  also 
is  termed  symmetric  (or  longitudinal)  motion,  and  the  associated  modes  of  vibration  are 
known  as  symmetric  modes.  In  flexural  motion,  the  average  displacement  is  transverse  to 
the  plane  of  the  plate.  This  is  also  termed  antisymmetric  motion,  and  the  associated  modes  of 
vibration  are  known  as  antisymmetric  modes.  As  a  result  of  the  combination  of  these 
modes  of  vibration,  any  disturbance  propagates  along  the  plate  at  a  speed  that  is  not  con¬ 
stant  but  depends  on  the  relative  ratio  of  wavelength  and  plate  thickness.  Knowledge  of  this 
dispersive  characteristic  of  wave  speed  in  plates  is  necessary  in  the  study  of  wave  scattering. 

The  exact  solutions  of  problems  of  wave  propagation  in  plates  or  layered  plates  of 
homogeneous,  isotropic,  and  elastic  materials  of  infinitely  extended  boundaries  have  been 
reported  in  the  literature.  The  analysis  is  very  complicated,  and  the  solutions  yield  relation¬ 
ships  between  frequency  and  wavelength  that  yield  an  infinite  number  of  dispersion  curves. 

If  the  plates  have  an  additional  constraint  such  as  a  stiffener  or  a  discontinuity  in  thickness 
or  in  material,  it  is  impossible  to  obtain  an  exact  solution.  As  a  result,  various  plate  theories 
have  been  established  as  means  of  approaching  such  problems. 

The  so-called  classical  theory  of  plate  vibration  is  based  on  the  Bernoulli-Euler 
theory  of  the  bending  of  a  beam  when  the  cross  section  of  the  beam  is  assumed  to  remain 
plane  and  perpendicular  to  the  neutral  axis  after  bending.  The  classical  theory  does  not 
describe  the  dispersion  of  the  first  antisymmetric  mode  of  vibration  and  hence  is  accurate 
only  for  a  thin  plate  at  low  frequencies  of  vibration.  The  Mindlin  theory  is  superior  when 
the  rotary  inertia  of  the  cross  section  and  the  shear  angle  of  the  neutral  plane  are  taken 
into  account.  In  fact,  the  importance  of  rotary  inertia  and  shear  angle  on  the  flexural 
motion  of  an  isotropic,  elastic,  homogeneous  beam  was  first  investigated  by  Timoshenko 
(ref  1 ).  Following  Timoshenko,  Mindlin  (ref  2)  formulated  an  approximate,  two-dimen¬ 
sional  theory  of  a  vibrating  plate  in  which  he  derived  a  single  equation  of  motion  involving 
transverse  displacement.  The  phase  velocity  of  straight-crested  waves  computed  from  his 
theory  was  seen  to  be  in  excellent  agreement  with  the  exact  theory  of  elasticity. 

For  a  two-layered  plate,  Jones  (ref  3)  used  the  exact  theory  to  compute  the  phase 
velocity  versus  the  thickness  of  the  two  layers.  Lai  (ref  4)  introduced  displacement  fields, 
which  include  the  shear  and  the  rotary  inertia  corrections  for  the  two  layers.  He  obtained 
an  analytical  solution  for  the  pressure  radiation  due  to  point  and  shear  force  excitations  on 
the  free  faces  of  the  plate  and  presented  the  numerical  values  of  this  solution. 


1  Timoshenko,  SP,  On  the  Correction  for  Shear  of  the  Differential  Equation  for  Transverse  Vibrations  of 
Prismatic  Bars,  Philosophical  Magazine,  vol  41,  p  744-746, 1921. 

2 

Mindlin,  RD,  Influence  of  Rotatory  Inertia  and  Shear  on  Flexural  Motions  of  Isotropic,  Elastic  Plates, 

J  of  Applied  Mech,  vol  18,  p  31-38, 1951. 

3  Jones,  JP,  Wave  Propagation  in  a  Two-Layered  Medium,  J  of  Applied  Mech,  vol  3 1 ,  p  2 1 3-222, 1 964. 

4  Lai,  JL,  Pressure  Radiations  from  an  Infinite  Two-Layered  Elastic  Plate  with  Point  and  Shear  Force 
Excitations,  J  of  Acoustic  Soc  of  America,  vol  53,  p  486-494, 1973. 
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Yu  (ref  5)  studied  the  sandwich  plate,  consisting  of  a  core  and  two  identical  faces. 

He  followed  the  analysis  of  Mindlin  (ref  2)  to  formulate  an  approximate  theory  that  includes 
corrections  for  shear  and  rotary  inertia.  In  the  limiting  case  where  either  the  core  approaches 
zero  thickness  or  the  two  layers  vanish,  the  sandwich  plate  becomes  a  single-layered  plate  and 
Yu's  results  become  identical  to  those  of  Mindlin's  work. 

In  many  applications  of  acoustics,  it  is  required  to  express  the  motion  of  a  two¬ 
layered  plate  in  a  single  equation  of  transverse  displacement,  but  no  such  equation  was  found 
in  the  literature.  This  study  is  aimed  at  correcting  that  deficiency.  The  present  analysis 
follows  very  closely  analyses  in  references  2.and  5.  Equations  of  motion  in  terms  of  transverse 
displacement  and  expressions  of  the  phase  velocity  of  straight-crested  waves  are  developed 
for  a  plate  consisting  of  two  elastic,  isotropic,  and  homogeneous  layers  of  arbitrary  thick¬ 
ness. 


FORMULATION 

Figure  1  shows  a  two-layered  plate  of  thickness  H.  The  two  elastic,  homogeneous, 
and  isotropic  layers  of  thickness  hj  and  hj  are  labeled  layer  1  and  layer  2,  each  with 
Young’s  modulus  Ej,  Poisson’s  ratio  vy  shear  modulus  #ij,  and  density  p-.  The  coordinate 
system  is  shown  with  plane  Oxy  coincident  with  the  interface  and  axis  Oz  normal  to  it.  The 
interface  contact  is  a  weld  and  the  layers  are  extended  infinitely  in  plane  Oxy.  To  simplify 
the  involved  algebra,  the  problem  is  restricted  to  the  plane  strain  condition ;  only  a  unit  length  of 
the  plate  needs  to  be  taken  in  the  y-direction.  The  displacement  in  the  y-direction  and  its 
derivatives  are  identically  zero.  The  free  faces  of  the  plate  are  shown  loaded  with  a  fluid 
that  exerts  normal  pressures  pzj  and  pz2  on  opposite  sides. 

Following  Mindlin  (ref  2)  we  first  define  the  equivalent  “plate  stress”  terms  and 
correlate  them  with  the  strain  components,  e^,  of  each  layer.  Then  we  assume  the  dis¬ 
placement  fields  ua  such  that  they  are  continuous  at  the  interface.  Finally,  we  use  the 
equations  of  motion  in  the  theory  of  elasticity: 

aa0,p  +  ~  • 

We  substitute  the  previously  assumed  displacement  fields  and  the  expressions  for  plate 
stresses  into  these  equations  of  motion  to  obtain  the  equations  of  motion  for  our  layered 
plate.  These  new  equations  of  motion  correlate  the  plate  stresses  and  the  assumed  displace¬ 
ment  fields.  By  using  the  plate  stress-strain  relationships,  we  obtain  a  set  of  equations  of 
motion  that  relate  the  assumed  displacement  fields  and  the  plate’s  elastic  properties. 


5  Yu,  YY,  A  New  Theory  of  Elastic  Sandwich  Plates-One  Dimensional  Case,  J  of  Applied  Mech,  vol  26, 
p  415-421,  1959. 
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FLUID 


PLATE  STRESS  COMPONENTS 

The  bending  moments,  in-plane  forces,  and  transverse  shearing  forces  are  defined 
respectively  as  follows: 


C  f 2 

MX1  J  °xxlzdz’  ^x2  “  J  °xx2z dz 

-hj  0 


Qxl  =  /  Tzxl  dz  >  Qx2  =  j1 


dz ,  Nx2  =  ^ 


Nxl  =  J  °xxl 
-hi 


rzx2  dz 


axx2  dz 


(1) 


where  the  subscripts  1  and  2  indicate  layer  1  and  layer  2. 
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From  Hooke’s  law: 

°a0  *  Xerr  6<*0  +  2tte*P  ’ 

We  can  write  the  stress-strain  relationship  as  follows: 


®xxi  "( 

)  exxi 

+  Xezzi 

°zzi  =  ( 

(  * 

|  czzi 

+  Xexxi 

Txzi  3 

^i^xzi 

and  Mi 

are  Lamfe’ 

s  constants  and 

=  _ 

Ev 

!  11  = 

E 

(1+v)  (l-2v)  2(l+»>) 

Eliminating  ezz  from  (2)  to  obtain  the  stress-strain  relationship: 


°xxi  "  j_  2  Cxxi  +  Xi+2Mi  °zzi 


Txzi  =  ^i^xzi 

Substitution  of  (3)  into  (1)  yields  the  following: 


>  i  Ei 


>i  xi 


d  ’  l  (W)'***  + 1 

d  =  J  **i^xzi 


°zziz  dz 


/i  Ei  ,bi  xi 

j^exxidz+J  ^7^°zzidz 


aj  *  Ki 


Since  there  is  no  concentrating  force  acting  on  the  plate  in  the  z-direction,  ozz j,  in 
general,  are  small.  Consequently  the  integrals  containing  a are  negligible  compared  to 
those  of  the  in-plane  stresses  oxxj.  It  is  therefore  assumed  that  these  integrals  are  equal  to 
zero,  which  leads  to  the  following  relations: 


8 


where  k  is  the  shear  correction,  introduced  exactly  as  in  Mindlin  (ref  2).  k  assumes  the  value 
of  jt^/12.  (Note  that  Mindlin  used  k?  instead  of  k.) 

PLATE  DISPLACEMENT  COMPONENTS 

Following  Mindlin  and  Medick  (ref  6)  and  Yu  (ref  5),  the  displacement  fields  can  be 
written  as  follows: 


Uj  (x,  z,  t)  =  u/0)(x,  t)  +  ZUj(1)(X,  t) 
Wj  (x,  z,  t)  =  w(x,  t) 


(6) 


where  the  superscripts  (0)  and  (1)  denote  the  order  of  displacement  and  the  subscript  i  de¬ 
notes,  as  before,  layer  1  and  layer  2.  u^^  represents  the  slope  of  the  cross  section  of  the 
layers  with  respect  to  the  z-direction.  In  writing  the  above  displacement  fields  we  made  two 
assumptions:  the  cross  section  of  each  layer  remains  plane  during  bending;  and  the  thick¬ 
ness  modes  (stretch  and  shear)  are  eliminated,  by  setting  wj  =  w  (x,  t).  Note  that  if  the  shear 
correction  were  neglected,  as  in  the  classical  plate  theory,  we  would  have 

Uj(1)(X,  t)  =  -  —  w(x,  t). 
ox 


The  displacement  fields  are  now  written  as  follows: 
u,  (x,  z,  t)  =  ^j(x,  t)  +  z\l>2(x,  t) 

u2(x,  z,  t)  =  \p  j (x,  t)  +  z\|/3(x,  t) 


(7) 


W|  =  W2(x,  z,  t)  =  w(x,  t) 

where  \p  j,  ^2,  ^3,  and  w  are  unknown  functions.  Notice  that  (7)  satisfies  the  condition  of 
continuity  of  displacement  at  the  interface  (z  =  0). 

Mindlin,  RD,  and  MA  Medick,  Extensional  Vibrations  of  Elastic  Plates,  J’of  Applied  Mech,  vol  26,  p  56 1  -569, 
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By  substituting  (7)  into  the  strain-displacement  relationship 

€ap  =  2  +  UP,^  ’ 

we  obtain 


3uj 

exxl 

3x 

11 

+ 

z^2 

3u2 

exx2 

3x 

11 

+ 

2  ^3 

3uj 

3wj 

'I'xzl 

3z 

II 

u 

+ 

1  ^2  +  w' 

3u2 

3w2 

?xz2 

3z 

II 

+ 

1  ^/3  +  w' 

(8) 


where  prime  ( ')  denotes  differentiation  with  respect  to  x.  Inserting  (8)  into  (5)  yields  the 
following: 


Mxl  = 


^l2 


hi3  \  2 


l-v{2 


2  *1  +  —  ^2  )  =  ~D1  *1  +  3  hl  D1  *2 


'-0  i  h22  h23  \  2 

Mx2=  - i  \—  ^\+-T  *3}  =D2*',  +-h2D2^ 

l-v« 


NX1  "  — 


hi2 


>-*V 


hi  ^l  — T~  ^2  )  =  TT  Di  ^l  ”  Di  ^2 


Nx2  = 


— £  \  h2  *1  +  ~  *3  j  =  °2  *1  +  D2  *3 


(9) 


Qxi  =  k  Ml  hj  (w'  +  \J/’2)  =  Gj  (w'  +  ^2) 

Qx2  =  K  ^2  h2  (w’  +  ^3)  =  g2  (w'  +  ^3) 
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where 


11 


Hence,  (11)  becomes 


M'xl  -Qxl  + 


[ZTxzl] 


0 

-h 


P\hl  v  P  lhl  v  A 

+  — —  \f/ j - —  *2  -  0  • 


1 


Similarly  for  layer  2,  we  have 


Mx2  *  Qx2  + 


P2h22  v 
—  *1" 


p2h23 

3 


^’3  =  0. 


02) 


(13) 


Note  that  the  last  terms  in  ( 1 2)  and  ( 1 3)  represent  the  rotary  inertia  of  the  cross  section  of 
each  layer. 

When  (10a)  is  again  integrated  over  the  thickness,  we  obtain,  for  layer  1, 


f  1  0  ••  Plhl2  •• 

N'xl  +  Txzl  -  P1M1  +— —  *2  -  0  . 
*■  *  -hj 


(14) 


and  for  layer  2, 


f  1  ^2  ••  ^2^2'  •• 

N'x2  +  Txz2j  -  P2h2*\ - -  *3  =  0  •  (15> 


(10b)  is  integrated  over  the  thickness,  yielding,  for  layer  1, 

Q'xl  +[  °zzl  ]  pjh,  w  =  0  ,  (16) 


and  for  layer  2, 


Q'x2+|°zz2  -P2h2w  =  0- 
’0 


(17) 
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The  boundary  conditions  of  this  problem  are  as  follows: 
At  z  =  -hj.  azz |  -  -Pzi  ,  rxzi  -  0  . 

At  z  ^2-  aa2  ~  ~  Pz2  >  Txz2  ~  ®  * 


At  z  -  0-  txzj  -  T"xz2  ’  °zzl  "  °zz2 
Applying  these  into  ( 1 2)  and  ( 1 3)  yields 

*>lhl2  v  *>lhl3 


M'xl  “  Qxl  +  2 

Mx2  "  Qx2  " 


- ^ —  ^2  =  0 


p2h22  ..  P2h23  •• 


*  1  ~ 


ip3  =  °  . 


(18) 


(19) 


We  now  apply  the  boundary  conditions  into  ( 1 4),  ( 1 5),  ( 1 6),  and  (17),  adding  ( 1 4) 
and  ( 1 5)  to  obtain 

P 1  h ! 2  ..  p2h22  .. 

N'x  -  (P  jh |  +  p 2^2^  ^1  +  — 2 —  ^2 - 2 —  ^3  =  ®  (20) 

and  adding  (16)  and  (17)  to  obtain 

Q'x  +  P  -  (P  jhj  +  p2h2)w  =  0  ,  (21) 

where 

Nx  =  Nxl+Nx2 
Qx  =  Qxl  +  Qx2 
p  =  Pzl  *  Pz2  • 

(18),  (19),  (20),  and  (21)  are  the  equations  of  motion  involving  the  plate  stresses  and  the 
assumed  displacement  fields. 

Substituting  the  expressions  for  plate  stress  components,  (9),  into  (18)  -  (21)  yields 
the  following: 

2  Plhl3  ••  Plhl2  •• 

-  Dj  4>"\  +  -  hjDj^'2  -  G|  (w'  +  \pj)  =  “  4*2 - 2 —  ^1  (22) 
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r 


D2  4>\ 


) 


2  P2^2  ••  P2^2 

+  "J  h2D2^3  ~  G2  (w’  +  ^3)  =  — 3 —  +  — 2 —  ^ 1 


i 


p2h22..  Pi^i2.. 

^1  ~  Dj  +  D2v!/3  =  — —  ^3 - 2 —  ^2 


i 


+  (p  1  h  1  +  P2h2)  $1 


(23) 


(24) 


(Gj  +  G2)  w"  +  Gj  \|/'2  +  G2$3  +  p  =  (p  1  h j  +  p2h2)  w  .  (25) 

(22)  -  (25)  are  the  equations  of  motion  connecting  the  assumed  displacement  fields  and  the 
elastic  properties  of  the  two  layers;  they  are  equivalent  to  equation  16  of  Mindlin  (ref  2) 
for  a  single-layered  plate.  In  practical  applications,  it  is  sometimes  required  that  the  equa¬ 
tions  of  motion  be  expressed  as  a  single  equation  of  w(x,  t).  This  can  be  done  by  eliminating 
\p  j,  ^2,  and  ^ 3  between  (22),  (23),  and  (24)  and  substituting  into  (25).  These  steps  are 
performed  next. 


>  1 

M 
*-  1 


RESULTS 

EQUATION  OF  TRANSVERSE  MOTION 

(22),  (23),  and  (24)  are  rewritten  as 

Lj  ^1  +  L2  ^2  =  L1 1  w 

L4*l  +l6  ^3  *  L12  w 
L7  i|/]  +  Lg  ^2  +  L9  1//3  =0 


where 


L,  =  -Dj  Vz  + 


Plhl  d]_ 
2  at2 


2 

3 


=  —  hjDj  V  -  Gj  - 


Plhl  32 
3  at2 


L A  =  Di  - 


p2h22  a2 

~  a? 


(22') 

(23') 

(24') 


I 
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M 

l7 


“  J  h2  d2  V2 


ii 

at2 


V2  -  (pjhj  +P2h2^ 


Lg  =-DjV2 


Ml  -  Gl* 

l12  =  g2V  . 


Hence  we  have 

A. 

[f]  *  Det  I  l4 

\h 


and 

r 

[£]  *3  =  Det  I  L4 

\l7 

where 


=  Det 


Miw 
M  2W 
0 

l2 

0 

l8 

M 

0 

Lg 


1 
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Substituting  the  above  into  (25)  yields 


+  £ 


12 


w 


(26) 


After  expanding  the  determinants  and  rearranging  the  terms  of  like  order,  we  obtain  the 
equation  of  motion  in  terms  of  w(x,  t): 

L,U  V2(h2L,C,  ♦  h^Gj-h^G,  -h,L4C2) 


where 


P  -  Pjhj  +P2h2 


and 


*  L1L4  [}  (h2Ll  '"I1'4)  +t(^  C2  +  h 

-j(G2L,2  +  G,L42,  +  2G,Gj^--j±) 


2 

1 


(27)  is  the  single  equation  of  motion  of  w(x,  t)  and  equivalent  to  equation  37  of  Mindlin 
(ref  2). 
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When  the  two  layers  are  identical, 


G1  -  G2 
Li  =  L4 
Pi  *  P2 

hi  =  h2 


After  some  algebraic  operations,  (27)  becomes 


f  /8  2  32  8  4  32 

[  1  V?'’1'’1  7?'Vh'Giv)  +TGi"i hi 

/  3°!  \  4hl 


w 


p  =  0  . 


(28) 


(28)  is  the  equation  of  motion  for  a  single-layered  plate;  it  becomes  identical  to  equation  37 
of  Mindlin  (ref  2)  if  h/2  is  substituted  for  hj.  We  substitute  the  expressions  of  Lj  and  L4  in¬ 
to  (27),  expand  it,  and  perform  some  algebraic  operations  to  yield  the  complete  form  of  the 
equation  of  motion: 


D  |  D2  (D  i  h2  +  D2hj)  (G  i  +  G2)  i 


+ 


-iGiGj(D,2  +  D22)-|D,D2-^U: 


Plhl\  P2h22hl\ 

- + - 1 


(Gi  +  G2) 


+  (P  i  h  j  +  P2h2)  (Djh2  +  D2hj) 


2  3 2  /  P 2h2" 

^lD'  — 


+  D- 


pihr 


(Djh2  +  D2hj)  (Gj  +  G2) 


"  D1D2^J  rI  r2  +  2^ 


G1G2 


,V4 


2_  <r_ 
9  at4 


DjD2  (pjhj  +P2h2) 


Pjhi  h2 


P2^2\\ 

+  “T7 
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+  (D1 
t(o 


P 2h22  plhl 

- +  D0 — — 


[\j  ^>lhl2h2  +  P2h22hl^ 

-> 


(Gj+Gj) 


Poho2  Pjhr  , 

— o —  +  d2 — o —  j  G°  lhl  +P2h2)(Dlh2  +  D2hl) 


p1h12p2h22 


(Djh2  +  D2hj)  (Gj  +  G2) 


32  f4 

+  — -  I  —  D j D2  (P 1 h  1  +P2h2)(rlG2  +  r2Gl) 

3t2  [3 

/  P2h22  ^  1 h  1 /  4  ,4  , 

+  ^1  2  +  D2  2  /  \Tri  +Tr2 +  2/  gig2 


+  -J  (Pjhj  +P2h2^  (Di2G2  +  D22  G j) 


+  y  G^VDj  p1h12  +  D2P2h22) 


+  V' 


IYd 

9  at6  IV 


p2h2 


P  Ihi 


1 — ~  +  D2 — ~  )(P\h\  +P2h2>\  2  + 


( 


p,hl2h2  P2b22hl 


PlP2hi2h22  /pjhj2h2  P2h22hA 

+ - ; -  \~~r-  +  2  /  (GI  + 


g2) 


PlP2h!2h22 


(P]h]  +  P2^2)  (Dj h2  +  ^2^1 )  I 


■M- 


-) 


p2b2z  p  ihr 

— “ — +  ^2 — 2 — /  (P 1  h j  +  p2h2>(riG2  +  r2Gl) 


PlP2hi2h22 


(j'l+j'2+2)C> 


+  — g!g2 


Plh, 


21  2 


P2h2 


G2 
2 1  2' 
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+  — -  (p |hj  +  p2h2)  l*1!2  +  GjD2P2^2  )j 


--f 

3t2  I 


2G1G2(p1h1  +P2h  2> 


\hl  h2  /  j 


8  A 

J_l_  (p1p2h12h22)(Pih1+P2h2>^lhrh2  +  p2h2  hl> 
36  gt8 


+  f— pjp2hj2h22  (Pjhj  +p2h2)  (rjG2  +  r2Gj) 

at6  l3 


+  “(Plhl  +P2h2) 


ill 

at4  i 


GjG2  (Pjhj  +P2h2) 


1 


+  G, 


P2h2 


21  2 


w  +  I  £  I  p  =  0  , 


with 


£1  =v6  |D1D2(h2D1+h1D2)j 


+  v 


1  a2  f  / p,hi2h2  P2h22hl\ 


(o, 


*>2h2' 


d2 


pihr 


(29) 


2  -^(h2DI  +h!D2)  -j  D,D2(r)G2  +  r2G1) 
J9S4  [/  P2h22  ,  Plhl2\/^ihl2h2  ^2h22h?\ 

+  V  |9^lVD,~  2~/V^  2  / 

+  ^>lP2h^l  .h_?Jj(h2Dl+hlD2)| 
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*7  (Dl<Vlhl2  ♦  D2G,/>2h22)j  +  2G,G2Ql  + 
+  j"  ^  PiP2hl2h22 <>lhl2h2 +  ^2h22hl>  JJ 

-  “4  |jPlP2hl2h22(r]G2  +  r2Gl> 

•h&y  4«.(¥)!i 


3^ 

-  2GjG2  (Pjhj  +P 2^*2^  — 2 


EQUATION  OF  PHASE  VELOCITY 

Consider  the  straight-crested  waves  propagating  in  the  x-direction.  The  fluid  loading 
term  is  not  needed  for  the  phase  velocity  equation.  We  assume  the  transverse  motion  of  the 
plate  as 


w  =  A0e‘<kx-kct)  . 


(30) 


After  substituting  (30)  into  (29)  we  expand  and  rearrange  terms  to  obtain  the  phase  velocity 
equation: 


£ 


8 


-^-T^tjRj  +  t2R2)  (R]3R22ti  +R]2R23t2> 


+  S 


-  3^T4(tlRl  +t2R2>  ^32  +  2/?22)Ri3R22tl  +  (P22  +  2p32)  Rj 

-  ^t4(‘1R1  +  Wl2>  [R13r22i1  *  +  Rl2R23t2«  I 

-  — T2k  (tjRj  +  t2R2)  |RiR2^R2rl^l2^2  +  Rlr2*l^ 

+  7'lRl^l2t7*2R23]j 
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+  *4  j^T4(t1R1+t2R2)[R13R22t1O324  +  202V> 

+  R12R23t2(034  +  2022032)] 

+  ~  T4  (tjRj  +  t2R2/3 , : 2)  [R,3R22t,  k  (fi3 2  +  2/322) 

+  R12R23t2»c(022  +  2032)] 

+  “T2  (tjRj  +  t2R2)  |RjR2  k  (^22  +  /332)  (R2rj/3j2t2  +  R^j) 

+  -y  ^22tlRl3#c^l2  +7^32t2R23Kj 
4T2,V[7'iV*7t22» 

+  K^^djRj  +t2R2)(t1R1  +t2R2) 

+  ^  j'^^(tlRl+t2R2)[Rl3R22tl%V  +  Rl2R2V2Vl 

-  “T4 (tjRj  +  t2R20(2)  [Rj3R22tj/c  (^24  +  2^22/332) 

+  R12R23t2K(^34  +  2022^32)] 

-  — T2k  (tjRj  +  t2R2)  |RjR2  (R2r  1  ^  1 2^2  +  RIr2l1^22^32 

+  -jR|V2V  +7R23<2%4] 

-7T2“2^2[f^2,'2R'4tf«32‘2V 

+  tf22  +«32)  t,t2R|2R22  <7  M  +  J  r2  +  2)1 


24  +  t,t2R,2R22(7r1  +7r2  +  2)] 
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-  K20l2(tlRl  +t2R2)0?22tlRl  +%2t2R2)[ 

+  jjgT4(t1R1+t2R2^12)[Ri3R22tl'‘^24^32)  +  Rl2R23t2K^22%4)| 

+  -^TV0I2(t12R14/324  +  t22R2V) 

+  ^-T2  t1t2R,2R22(c2/312/322/332  (  ~t1  +yr2  +  2)J  =0  ,  (31) 

where 

f  =  c2/cs,2  ,  p2  =  cs22/cs]2  , 

022  =  cpi2/csi2  ,  032  =  cp22/cs]2 

hj  _  h2 

Fl  h2  ’  F2  hj 

hl  h2 

h|  +  h2  =  H  ,  Rj  ~  -  r2  ~  H 

T  =  kH, 

Pi  P  2 

ti  =  -  ,  t2 - - - 

1  P1+P2  P1+P2 


and 

c  =  cj/k  =  phase  velocity  of  the  straight-crested  waves 

csi »  cs2  =  shear  wave  velocity  csj  =  (Mj/Pj)1^2  in  layer  1  and  layer  2 

Cpl ,  cp2  =  longitudinal  plate  wave  velocity  cpj  =  [  Ej/pj  (1  -  i»j2)  ]  ^  *n  layer  1  and  layer  2. 

(31)  is  the  equation  of  the  phase  velocity.  Examination  shows  four  branches  of  the  phase 
velocity  curves  versus  T,  the  wavenumber-thickness  product.  It  is  readily  verified  that  when 
the  two  layers  are  identical,  (3 1 )  becomes  the  equation  for  a  single  plate  of  the  Timoshenko- 
Mindlin  model. 

NUMERICAL  RESULTS 

The  numerical  results  of  phase  velocity  curves  were  computed  for  a  component  plate 
of  a  steel  layer  (layer  1 )  and  a  copper  layer  (layer  2).  The  elastic  properties  of  these  layers 
are  as  follows: 
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Steel 


Ej  *  2.168663  X1011  N/m2 
vi  ■  0.283629 
pj  =  7.8  X103  kg/m3 
Cpj  =  5498.7  m/s 
csj  =  3290.9  m/s 
Cri  =  3044.9  m/s 
Copper 

E2  =  1.22  X1011  N/m2 
(>2  =  0.33 

p2  =  8.9  XI 03  kg/m3 
Cp2  =  3922.1  m/s 
cs2  =  2270.1  m/s 
cr2  =  21 15.8  m/s 

crj  are  the  Rayleigh  wave  velocities. 

The  normalized  phase  velocity  £  =  c/csj  is  plotted  versus  the  wavenumber-thickness 
product,  T  =  kH.  Appendix  A  shows  the  resulting  curves  for  nine  thickness  ratios  Rj  =  h  j/H. 
The  following  comments  pertain  to  those  curves. 

As  T  increases  from  zero,  curves  1  (first  antisymmetric  mode  of  motion  )  increase 
rapidly  until  T  *  2  and  gradually  reach  an  asymptotic  value  lying  within  the  interval 
cR2/,csl  “  CR l / cs  1  an^  increasing  with  thickness  ratio  R  j.  This  can  be  explained  easily  on 
the  following  physical  ground.  For  the  limiting  case  of  a  single  layer,  curves  1  attain  the 
Rayleigh  wave  speed  for  large  T;  for  the  case  of  two  layers,  the  asymptotic  value  of  curves  1 
tend  to  reach  the  Rayleigh  wave  speed  of  the  thicker  layer.  Furthermore,  the  shape  of 
curves  1  is  very  similar  to  that  of  the  exact  solution  (ref  3). 

Curves  2,  representing  the  phase  velocity  of  the  second  antisymmetric  mode  of  motion, 
reach  asymptotically  the  plate  wave  speed  of  the  soft  layer,  in  this  case  1.19.  When  T  =  0, 
the  starting  points  of  these  curves  increase  from  Cp2/csj  to  Cpj/csj,  or  from  the  plate  wave 
speed  in  copper  to  that  in  steel  as  Rj  increases  from  0.1  to  0.9.  Similarly,  curves  4,  represent¬ 
ing  the  phase  velocity  of  the  fourth  antisymmetric  mode  of  motion,  reach  asymptotically  the 
plate  wave  speed  of  the  hard  layer  -  in  this  case  steel,  with  Cpj/cs|  =  1.67.  The  asymptotic 
values  of  curves  2  and  4  do  not  change  with  varying  R  j. 

Curves  3,  representing  the  phase  velocity  of  the  third  antisymmetric  mode  of  motion, 
are  seen  to  wander  within  the  band  delimited  by  curves  2  and  4.  The  asymptotic  values  of 
curves  3  behave  similarly  to  those  of  curves  1,  shifting  toward  the  plate  wave  velocity  of  the 
thicker  layer  as  Rj  increases  from  0.1  to  0.9. 
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CONCLUSIONS 


An  approximate  theory  of  motion  of  a  two-layered  plate  of  elastic,  isotropic,  and 
homogeneous  materials  was  developed.  Shear  and  rotary  inertia  corrections  were  included 
in  the  theory.  Equations  of  motion  and  of  the  phase  velocities  of  straight-crested  waves 
were  derived.  In  the  limiting  case,  these  equations  were  shown  to  be  identical  to  those  of 
single-layered  plate.  Numerical  results  showed  that  dispersive  curves  of  the  phase  velocities 
versus  the  wavenumber-thickness  product  of  the  first  four  antisymmetric  modes  of  motion 
behave  very  similarly  to  those  of  the  exact  solution  developed  by  Jones  (ref  3).  The  phase 
velocity  of  the  lowest  mode  reaches  an  asymptotic  value  lying  within  the  individual  Rayleigh 
wave  speeds  of  the  two  layers.  The  asymptotic  values  of  the  phase  velocities  of  the  second 
and  the  fourth  modes  are  recognized  as  the  plate  wave  speeds  of  the  two  layers.  The  curve 
representing  the  third  mode  is  seen  to  wander  within  those  asymptotic  values  of  the  second 
and  fourth  modes. 


RECOMMENDATIONS 

1.  Employ  this  model  to  investigate  the  behavior  of  the  scattered  signal  (target 
strength)  from  a  two-layered  plate.  A  useful  approach  is  to  assume  the  plate  to  be  loaded 
with  water  on  one  side  and  air  or  water  on  the  other  side,  then  to  compute  and  plot  the 
transmission  and  reflection  coefficients  of  an  incident  acoustic  wave.  The  procedure  for 
performing  the  calculations  is  very  similar  to  that  of  Graff  et  al  (ref  7).  The  variation  of 
the  coefficients  versus  the  angle  of  incidence  and  the  wavenumber-thickness  product  yields 
information  about  target  strength  and  its  variation  at  various  modes  of  plate  vibration. 
Furthermore,  since  the  coefficients  can  be  computed  from  the  exact  theory,  a  quantitative 
comparison  can  be  made  and  the  validity  of  this  approximate  model  can  be  examined. 

2.  Apply  this  model  to  a  composite  plate  of  one  elastic  layer  (base  plate)  and  one 
viscoelastic  layer  (coating  layer),  to  calculate  target  strength.  Impedance  discontinuities 
along  the  plate  such  as  stiffening  ribs  and  varying  thickness  can  be  handled  by  appropriately 
modifying  the  model. 

3.  Note  that  this  model  can  also  be  used  where  the  viscoelastic  layer  contains 
many  inclusions  and/or  cavities,  provided  the  “equivalent  effective  elastic  properties”  of 
the  viscoelastic  layer  are  calculated  and  used. 

4.  If  experience  with  a  single-layered  plate  is  any  indicator,  this  model  should  be 
improved  to  provide  better  resolution  of  the  scattered  signal  of  acoustic  waves  at  grazing 
incident  angles,  possibly  by  combining  it  with  the  “Lyamshev  Theory”  on  symmetric  modes 
of  vibration. 


y 
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APPENDIX  A:  PLOTS  OF  NORMALIZED  PHASE  VELOCITY 
VS  WAVENUMBER-THICKNESS  PRODUCT 


Each  figure  consists  of  four  dispersive  curves  of  dimensionless  phase  velocity  versus 
the  wavenumber-thickness  product,  for  a  given  thickness  ratio  Rj  =  h  j/H.  The  four  curves 
are  labeled  1-4  to  represent  the  first  through  the  fourth  antisymmetric  mode  of  motion. 
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Figure  A7.  Rj  =  0.7 
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